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Dephasing of qubits by transverse low-frequency noise 
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We analyze the dissipative dynamics of a two- level quantum system subject to low-frequency, e.g. 1/f 
noise, motivated by recent experiments with superconducting quantum circuits. We show that the effect of 
transverse linear coupling of the system to low-frequency noise is equivalent to that of quadratic longitudinal 
coupling. We further find the decay law of quantum coherent oscillations under the influence of both low- and 
high-frequency fluctuations, in particular, for the case of comparable rates of relaxation and pure dephasing. 
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Recent experiments with superconducting 
Josephson-junction circuits ^ |2] demonstrated 
quantum coherent oscillations with a long decay time 
and a quality factor up to ~ 10 4 . These experiments, 
on one hand, probe coherent properties of Josephson 
qubits (quantum bits) and demonstrate their potential 
Tor applications in quantum computing and quantum 
communication. On the other hand, they may be 
viewed as a probe of the noise mechanisms in the 
devices studied. 

For the description of the dynamics of a two-level 
system (qubit, spin) subject to weak short-correlated 
noise one may use the Bloch equation, known from the 
NMR studies, which describes exponential relaxation of 
the longitudinal spin component and dephasing of the 
precessing transverse spin component (here and below 
we use the spin-1/2 language to discuss the dynamics). 
This description is valid as long as the correlation time 
of the noise is short compared to the typical dissipative 
times T\, Tg. However, in Josephson-junction qubits the 
low-frequency noise is strong. These low-frequency fluc- 
tuations are correlated over distant times, and special 
treatment of their influence on a qubit is needed. They 
could lead to complicated decay laws [31 EI El IH]- In 
Ref. [7] the influence of low-frequency fluctuations non- 
linearly coupled to a qubit was analyzed; this analysis is 
relevant for operation at the so-called optimal operation 
points pp. Here we extend this analysis to account for 
the effect of transverse fluctuations also present at opti- 
mal points. While our discussion applies to an arbitrary 
dissipative two-level system, for illustration we consider 
the Josephson charge qubit, similar to that studied in 
the experiment . We begin by discussing this system 
and the relevant noise sources and then proceed to the 
analysis of dephasing in general and at optimal points. 



n 



C V„ 



FIGQ The simplest Josephson charge qubit 

Dissipative dynamics of a Josephson charge 
qubit. The simplest Josephson charge qubit is the 
Cooper-pair box shown in Fig^ [SI- It consists of 
a superconducting island connected by a dc-SQUID 
(effectively, a Josephson junction with the coupling 
_Ej(<I> x ) = 2E® cos(7r$ x /$o) tunable via the magnetic 
flux <£> x ; here <£> = hc/2e) to a superconducting lead 
and biased by a gate voltage V g via a gate capacitor C g . 
The Josephson energy of the junctions in the SQUID 
loop is Ej, and their capacitance Cj sets the charging- 
energy scale E c = e 2 /2(C g + Cj), Cj = 2C?. At low 
enough temperatures single-electron tunneling is sup- 
pressed and only even-parity states are involved. Here 
we consider low-capacitance junctions with high charg- 
ing energy Ec 3> Ej. Then the number n of Cooper 
pairs on the island (relative to a neutral state) is a 
good quantum number; at certain values of the bias 
V g w Vdog = (2n+ l)e/C two lowest charge states n and 
n + 1 are near-degenerate, and even a weak Ej mixes 
them strongly. At low temperatures and operation fre- 
quencies higher charge states do not play a role. The 
Hamiltonian reduces to a two-state model, 
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in the basis J.) = |n), ||) = |n+l); here E c h(V g ) = 
2e(Vg - Vdcg)C g /(Cj + C g ). The Hamiltonian {TJ can 
be controlled via the gate voltage V g and the applied 
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flux <I> X , this allows one to manipulate the qubit's state 
and perform quantum logic operations. To read out the 
final quantum state one has to couple the qubit to a 
quantum detector, e.g. a single-electron transistor [S]. 

Quantum bits are inevitably coupled to fluctuations 
in the environment (bath) . This destroys the coherence 
of the qubits' dynamics. To slow down the dephasing 
the coupling should be made as weak as possible. In 
solid-state systems decoherence is potentially strong due 
to numerous microscopic modes. In Josephson qubits 
the noise is dominated by material-dependent sources, 
such as background-charge fluctuations or variations of 
magnetic fields and critical currents, with power spec- 
trum peaked at low frequencies, often 1//. A further 
relevant contribution is the electromagnetic noise of the 
control circuit, typically Ohmic at low frequencies. The 
1// noise appears difficult to suppress and, since the de- 
phasing is dominated by low-frequency noise, it is par- 
ticularly destructive. On the other hand, Vion et al. 
showed that the effect of this noise can be substantially 
reduced by tuning the linear longitudinal qubit-noise 
coupling to zero (in a modified design; they also sup- 
pressed the coupling to the quantum detector to mini- 
mize its effect on the qubit before the read-out). This 
increased the coherence time by 2-3 orders of magnitude 
compared to earlier experiments. 

Of special interest is the analysis of the slow dephas- 
ing at such an optimal point. On one hand, comparison 
of theory and experiment may verify our understand- 
ing of the physics of the device studied as a dissipative 
two-level system. Further, from the analysis of the de- 
phasing time scale and the decay law one may extract 
additional information about the statistical properties 
of the noise. On the other hand, understanding of the 
dissipative processes should allow their further suppres- 
sion in future qubit designs. 

Part of the noise (inch the background-charge fluc- 
tuations) can be thought of as fluctuations of the gate 
voltage and another part as fluctuations of the control 
flux $ x . It is convenient to discuss the effect of, e.g., the 
voltage noise V g = V® + Y (t) in the qubit's eigenbasis: 

H = --[AE o z + (Y(t)(- sin 77 & x + cost? & z )} , (2) 

where the level splitting AE = (^(Vg ) 2 + E 2 ) 1 / 2 and 
the angle between the static and fluctuating 'magnetic' 
fields is given by tan?7 = Ej/E c h(Vg). We expanded 
the variation of E c ^ in Y to the linear order. Consider 
first the effect of weak short-correlated noise (with cor- 
relation time shorter than the dissipative times; this 
includes the finite-temperature Ohmic noise). In this 
case one can use the lowest-order perturbation theory 



and finds that the spin dynamics is described by the 
Bloch equations, known from NMR. The interlevel tran- 
sitions are induced by the transverse fluctuations oc sin 77 
and give the relaxation time 1 jT\ = £ 2 sin 2 r/Sy (^ = 
AE)/2; the dephasing time is 1/T 2 = 1/(2T X ) + 1/T 2 *, 
where the pure dephasing is induced by the longitudinal 
noise cx cos ij and gives l/X^* = ( 2 cos 2 nSy (w = 0)/2 
(here the noise power S Y (t) = (l/2)([Y(t),Y(0)]+); we 
set h = 1). The effect of the magnetic-flux noise can 
be analyzed similarly. For Josephson qubits these ex- 
pressions give good estimates for the measured relax- 
ation times but do not suffice to describe the dephas- 
ing. Indeed, the expression for T 2 * cannot be used for 
strong longitudinal low- frequency, e.g. 1// noise; still 
it indicates that dephasing is strong. In first experi- 
ments E] dephasing times in the range of fractions to 
a few nanoseconds were achieved. Tuning to an optimal 
point extended the coherence time to ~ 1 /zs . 
Dephasing at optimal point. We illustrate our dis- 
cussion of decoherence at an optimal point by consid- 
ering a qubit deep in the charge limit, although in the 
device of Ref. Eq and Ej were comparable (in which 
case two lowest eigenstates, which form the qubit, are 
no longer charge states). Using two control parameters 
V g and $ x one can tune the longitudinal linear couplings 
to the charge and flux noise to zero: For instance, for 
the system tuning the gate voltage to the degener- 
acy point E c h(Vg) = yields cos?] = 0. Further, tuning 
<f> x to the point of maximal i?j(<I> x ) also suppresses the 
linear coupling to the flux fluctuations $ x = $ x + X(t). 
Thus, at this optimal point the Hamiltonian reads: 

H=-~[AEa z + XX 2 a z +(Ya x ], (3) 

where we left only the leading fluctuating terms. 

The quadratic londitudinal low-frequency noise 
XX 2 a z may result in an unusual dephasing law (with a 
power law crossing over to exponential decay) due to 
strong higher-order contributions (7j- Here we discuss 
the effect of the transverse noise CYa x . It can lead to 
relaxation processes and contribute to pure dephasing 
in higher orders. Thus in the analysis of dephasing one 
needs to account for both XX 2 and C,Y terms. 

The effect of the low- frequency (w <C AE) transverse 
noise can be treated in the adiabatic approximation: 
we diagonalize © to -a zy /(AE + XX 2 ) 2 + ((Y) 2 /2 w 
-& z [AE + XX 2 +( 2 Y 2 /(2AE)}/2, thus the low-w trans- 
verse noise contributes to pure dephasing. In general 
higher-frequency fluctuations are also present and in- 
duce relaxation. If the relaxation is much slower than 
the pure dephasing, one may neglect its contribution to 
the total dephasing. If the relaxation is much faster, 
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it dominates the decoherence; in this limit its rate 
£ 2 Sy(AE)/2 is given by the golden rule. However, of 
special experimental interest is the situation with 
comparable relaxation and pure-dephasing time scales. 
We analyze whether evaluation of each of them is influ- 
enced by the other in this case, that is whether the low- 
and high-w contributions interfere. In particular, we 
expect strong higher-order contributions to the pure 
dephasing due to strong low-w noise. Does it also con- 
tribute to relaxation? In the lowest order the relaxation 
is due to transitions with emission of a single resonant 
bath excitation; can instead a near-resonant excitation 
be emitted accompanied by low-frequency excitations? 
Here we show how the dephasing and relaxation laws 
and time scales can be obtained. 

Dephasing by transverse noise. We begin with a 
discussion of purely transverse noise, A = 0. We focus 
on the long-correlated noise (slow decay of (Y(0)Y(t))), 
i.e. on the noise power peaked at low and smooth at 
high frequencies. In our analysis below we assume weak 
dissipation: pure dephasing and relaxation slower than 
the osillations, T <C AE, where T represents the total- 
dephasing time scale. This limit is of primary interest 
for the circuits that realize qubits. Further, below for 
illustration we consider a source of gaussian noise Y(t), 
which can be characterized by its second correlator, but 
our major conclusions persist in more general situations. 

Our discussion is based on the analysis of the evolu- 
tion operator of the qubit dynamics using the 'real-time' 
Keldysh diagrammatic expansion in the qubit-bath cou- 
pling (this approach |10j is useful, since the spin degree 
of freedom does not satisfy the Wick theorem; it reminds 
the approach of Ref. [H]) . We begin by showing that the 
subleading-order effects of the low-frequency transverse 
noise reduce to the lowest-order contribution of longitu- 
dinal quadratic noise (this can also be seen from the adi- 
abatic approximation but our derivation indicates the 
diagrams, important in the discussion below). 

In the diagrams the horizontal direction explicitly 
represents the time axis. The solid lines describe the 
unperturbed (here, coherent) evolution of the qubit's 
2x2 density matrix p, exp(— iLot)9(t), where Lq is the 
bare Liouville operator (this translates to l/(u> — iLq) in 
the frequency domain). The vertices are explicitly time- 
ordered; each of them contributes the term QY<j x t z /2, 
with the bath operator Y(t) and the Keldysh matrix 
t z = ±1 for vertices on the upper/lower time branch. 
Averaging over the fluctuations should be performed; 
for gaussian correlations it pairs the vertices as indicated 
by dashed lines in FigEl each of the lines correspond- 
ing to a correlator (YY). Fig|3 shows contributions to 

(2) 

the second-order self-energy H-j.^^ (here ij =H label 
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FIG. [5] Second-order contributions to the self-energy 
£fj<_1-j_. Other terms are obtained by shifting an even 
number of the vertices in a or b to the lower branch. 

four entries of the qubit's density matrix). The term in 
FigE^L gives 



|) ldndr 2 {Y(t)Y(t')){Y(n)Y(r 2 ))e iAE ^-^ , (4) 



with integration over the domain if < T2 < t\ < t. 

After the summation over vertex positions on the 
lower/upper branches in Figs|2tjb, we evaluate the be- 
havior of the Fourier-transformed self-energy in the 
vicinity of the level splitting, at u> — —AE — u>' + iO, 
where ui' <C AE. If the integral is dominated by low 
frequencies, we find: 

x [<^ 2 w><n> + <y?- u >)<y?u)] = -\syM) , (5) 



where Sy 2 is the noise of 

Y 2 = ( 2 Y 2 /(2AE) . 



(6) 



The result (JSJ) coincides, as expected, with the lowest- 
order contribution of the term —Y2 & z /2; note that the 
left and right vertex pairs in Fig|2lcan be viewed as com- 
posite vertices corresponding to —Y2 <J z /2. Below we 
demonstrate that similar reduction occurs in every or- 
der of the perturbative expansion. Specifically, we show 
that the decay of the off-diagonal entry of the density 
matrix is /o-fi (*) = PuW 'exp(— t/2T\) -7 v (i), where the 
relaxation time is given by the golden rule and the pure 
dephasing term j v (t) is the same as for the longitudinal 
fluctuations — Y2 &z/2 (analyzed in Ref. 0). 

To demonstrate this we consider the diagrammatic 
calculation of the evolution operator for the density ma- 
trix. We begin by evaluating the evolution of the off- 
diagonal entry p, (the phase dynamics and dephasing), 
and then discuss relaxation (evolution of pj, pj). 

The expansion of the propagator contains vertices 
on the horizontal solid lines, representing the Keldysh 
contour. For a given number and ordering of the ver- 
tices, one is to integrate over their time positions, and 
then add up all diagrams. Consider the dependence of 
the integrand on the time position of a vertex. This de- 
pendence includes fast oscillations with frequency ±AE, 
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since the vertex flips the spin and changes the energy 
of the bare hamiltonian — AE a z /2, and a much slower 
dependence of the dashed line. Thus the integrand is 
a fast oscillating function of the time positions of each 
vertex. Since integration is typically performed over 
time range much wider than the period of oscillations 
(at all times of interest for the analysis of dephasing; 
this range is ~ 1/r at t ~ 1/r), the contribution of 
the most part of the integration space (with vertices' 
time positions as coordinates) is strongly suppressed by 
fast oscillations. However, in certain directions in this 
space, in which pairs of vertices with opposite oscilla- 
tion frequencies ±AE move together, the variation is 
slow, and the respective domains dominate the integral. 
One can arrive at this conclusion, and determine the 
dominant domains, by considering the evaluation of a 
particular diagram: the integral over the time t of a 
vertex J g(t) exp(iAEt)dt is taken between the posi- 
tions a, b of the neighboring vertices. Since g{t) is slow 
on scale 1 / AE, the oscillatory integral is dominated by 
the boundary terms, g(t) exp(iAEt)\ b a / (iAE) . One can 
say that the vertex t joins one of its neighbors, and later 
one integrates over the vertex-pair position, a or b. 

One can continue this process, integrating at each 
step over time positions of unpaired vertices (or clusters 
with an odd number of vertices and hence oscillatory be- 
havior), if any are still present. Finally one arrives at a 
situation where all vertices are paired, and the depen- 
dence of the integrand on pair time positions is slow (the 
exponentials exp(±«A_Ei) for two paired vertices com- 
pensate each other). The integral in each domain with a 
fixed time ordering is dominated by the boundary terms, 
that is the terms with paired vertices. Thus we elimi- 
nate the high-frequency (~ AE) behavior, and now can 
evaluate the propagators using the diagram technique 
with ingredients that are slow (without oscillatory de- 
pendence on their time position): 'double' vertices with 
two dashed tails in FigJSK, and dashed lines connect- 
ing these tails (cf. the examples in FigEb). Although 
4-, 6- and further 2n-fold clusters also form slow ob- 
jects, their creation requires additional constraints on 
the vertex times (compared to building n pairs) and the 
respective integration domain is much smaller; thus the 
contribution of such clusters is of higher order in Y / AE. 

A closer inspection of the spin dependence and the 
Keldysh two-branch structure reveals that in each pair 
both vertices are located on the upper or both on the 
lower branch (for vertices located on different branches, 
two terms with different time orderings cancel each 
other unless the vertices are linked by a dashed line; 
such a term appears in the analysis of relaxation but 




FIG. El a. Double vertices with \aw-u tails, which ap- 
pear in the evaluation of dephasing. b. Examples of 
clusters built out of them 0. c. A low-w object with a 
high-frequency dashed line. The relaxation process in 
e also contributes to dephasing as shown in d. 



not of the dephasing), and they indeed effectively cor- 
respond to the term —Y% a z /2 in the Hamiltonian. 

So far we constructed slow composite objects pay- 
ing attention only to the oscillations of the solid lines 
in the diagrams and assuming very slow dashed lines, 
i.e., neglected the higher-frequency noise. In fact, one 
can construct another slow object shown in FigEh, if 
the respective oscillations of the solid lines are com- 
pensated by the dashed line from this vertex. In other 
words, in the frequency domain, one constrains the fre- 
quency of the dashed line to be AE (or — AE, depepnd- 
ing on the direction of the spin flip at the vertex) . The 
dashed lines from such objects pair up, and the integral 
w.r.t. their relative position is dominated by small sep- 
arations, St ~ l/AE. Thus one finds the slow object 
of Fig[3Ji, two vertices linked by a dashed line at fre- 
quency AE; it describes the relaxational contribution 
to dephasing exp(— t/2T\), where 

l/T 1 = ( 2 S Y (AE)/2. (7) 

In similar clusters of higher order additional constraints 
strongly limit the integration domain. Note that the 
object in FigOJi involves the weak noise at a high fre- 
quency AE, unlike those in Fig|3^, but it is still relevant 
since the lowest-order term in the upper part of Fig[3|3 
is imaginary and does not contribute to dephasing. 

Similarly, we analyze the relaxation of the diagonal 
entries p j , pj. The new slow ingredient in this analysis is 
shown in Fig|3^. As for the composite objects in Fig|3Ji, 
the terms with these objects located on the upper and 
lower branches cancel each other, due to different signs 
ascribed to them in the Keldysh formalism (in contrast, 
in the analysis of the evolution of pj displacing a vertex 
from one branch to the other flips the spin thus yield- 
ing an additional sign change rather than cancellation) . 
Hence the relaxation is given by the terms in Fig[3|i,e 
and Eq. Q . We find that the strong low- frequency noise 
does not influence the relaxation rate. 
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Discussion. We focused on the effect of the purely 
transverse noise. One can verify that in the presence of 
the longidsutinal fluctuations —XX 2 a z /2 the reduction 
persists: the relaxation is still given by Eq. Q, and the 
dephasing can be found by considering the longitudi- 
nal noise - [XX 2 + ( 2 Y 2 /(2AE)]a z /2. For uncorrected 
fluctuations X{t) and Y(t) their effects just add up (this 
would happen in the charge limit Eq 3> Ej for the qubit 
in FigQ]at the degeneracy point, where the charge noise 
is transverse, and at the proper flux bias, where the flux 
noise is quadratic longitudinal). In the experiment 
Eq and Ej were comparable, hence both charge and flux 
noise contributed to longitudinal and transverse fluctua- 
tions making them correlated; this should be taken into 
account but does not complicate the analysis. 

We considered slow fluctuations. For short- 
correlated noise the double vertices of Fig|3Ji do not 
contribute, and one recovers the Bloch equations. 

Furthermore, we illustrated our analysis by an ex- 
ample of gaussian noise. Such fluctuations are in- 
deed encountered: The low-frequency noise, e.g. the 
background-charge fluctuations in Josephson circuits, is 
possibly produced by a collection of bistable fluctuators 
(or discrete system with more states). With a proper 
wide distribution of their parameters (couplings to the 
qubit, switching rates) they produce a smooth 1// noise 
power. If the qubit is affected by many fluctuators, with 
a dense distribution in the parameter space, due to the 
central limit theorem one expects gaussian noise. In 
some Josephson devices sharp noise features indicate 
that a few fluctuators dominate and the resulting noise 
is non-gaussian (dephasing by such fluctuators was stud- 
ied, e.g., in Ref. [Hj)- We emphasize that our reduction 
applies also to the analysis of these systems. Indeed, 
our derivation used only the fact that the fluctuations 
are slow. Thus one can still build the diagrams from 
the slow objects constructed above but for non-gaussian 
fluctuations the dashed tails of the vertices may join not 
in pairs but also in larger bunches. After the reduction 
to quadratic longitudinal noise one may use other, non- 
diagrammatic, ways to analyze its effect. 

Our results are relevant for the analysis of the ex- 
periment of Ref. PP, in which the measured relaxation 
and dephasing times were comparable. The prediction 
of a specific decay law requires a detailed knowledge of 
the noise power of charge and flux fluctuations. It can 
be acquired via measurements away from the optimal 
point, as indicated in Ref. p. 
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